First Constraints on Source Counts at 350 Microns 
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We have imaged a ~6arcminute 2 region in the Bootes Deep Field using 
the 350 /xm-optimised second generation Submillimeter High Angular Resolution 
Camera (SHARC II), achieving a peak la sensitivity of ~5mJy. We detect three 
sources above 3a, and determine a spurious source detection rate of 1.09 in our 
maps. In the absence of 5a detections, we rely on deep 24 /im and 20 cm imaging 
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to deduce which sources are most likely to be genuine, giving two real sources. 
From this we derive an integral source count of 0.84^ '^ sources arcmin -2 at 
S > 13mJy, which is consistent with 350 /im source count models that have an 
IR-luminous galaxy population evolving with redshift. We use these constraints 
to consider the future for ground-based short-submillimetre surveys. 

Subject headings: infrared: galaxies - submillimeter: galaxies - galaxies: star- 
burst - galaxies: high-redshift 



1. INTRODUCTION 

Towards the end of the last decade, a new population of submillimeter-selected galaxies 



at the James Clerk Maxwell Teles cope (jSmail. Ivison &: Blainl 119971 ; iHughes et al. 



Holland et al. 


1999) 


iHuehes et al. 


1998 



Barger et al.lll998l ; lEales et al.lll9991 ). Extensive follow-up obser yations constrained the pop 



ulation to be mainly m assive star forming gal axies (see, e.g., 



Chapman et al. 



bx et al. 2002: Blain et al. 



20051 ). with the bulk of the 



20021 : iBorvs et all 120031 ) at high redshift (z~2; 
luminosity emitted in the restframe far-IR — although the detection of these sources had 
not been predicted by semi-analytical hiararchical models (fo r instance, contrast the order 
of magnitude s pread in the model s of iGuiderdoni et al.lll997l with the SCUBA-constrained 
single model of iGuiderdonil Il998af ) . SMGs can be considered the high redshift counterparts 
of the luminous and ultraluminous infrared galaxies (LIRGs and ULIRGg^) foun d in the local 



universe (the majority select e d with the Infrare d Astronomical Satellite, IRAS; ISoifer et al. 
19841 : 1 Joseph fc WrkStlll985l : ISoifer et al.lll987h . 



The first surveys with SCUBA paved the way for m any similar suryeys using other sub 



millim etre detectors (e.g., MAMBO, BOLOCAM, etc., iBertoldi et all l2000l : lLaurent et al. 



20051 ). These surveys were mostly limited to long submillimetre through millimetre wave- 
lengths (500-1300 /im), tracing emission on the long wavelength side of the peak at typical 
redshifts. The commissioning of the 350 /xm-o ptimised Second G eneration Submillimeter 
High Angular Resolution Camera (SHARC II; iDowell et al.l 120031 ) at the Caltech Submil- 
lime ter Observatory (CS O), currently the largest ground-based submillimetre bolometer ar- 
ray (IMoseley et al.l 12004 ). provided a feasible opportunity to carry out a blind survey in this 
waveband. 



defined by 10 11 L Q < L m (8-1000 /im) < 10 12 L Q and L IR (8-1000 /urn) > 1O 12 L respectively 
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Using SHARC II, we targeted a ~ 6 arcminute 2 region in the Bootes Deep Field (Ide Vries et al 



20021 ) for a blank, deep survey, which was designed to select LIRGs and ULIRGs through 
their far-IR thermal dust emission (measured near the peak of the spectral energy distribu- 
tion (SEP)) at 1 < z < 3, the epoch of peak cosmic star formation rate density (see, e.g., 
Hopkins fc Beacomll2006l ). The survey, achieving a peak la sensitivity of ~5mJy, produced 
a promising result as we reported the discovery of the first galaxy selecte d at 350 jum (SMM 
J143206.65+341613.4, also named Short Submillimetre Galaxy 1 (SSG l): lKhan et alJl2005h . 



The discovery of SSG 1 raises a number of questions regarding the nature of galaxies 
detected in the short-submillimetre wavebands (200-500 yum). Given the demanding obser- 
vational requirements (good 35 nm atmospheric tr ansmission from Mauna Kea is ~30% 
as opposed to ~80% at 850 /zm; ISerabyn et al.lll998l ). short-submillimetre surveys would be 
a poor use of ground-based telescope time if they traced the same population as did long- 
su bmillimet re s urvey s . However, follow-up 1.2 mm imaging appears to confirm the assertions 
of iKhan et al.l (120051 ). that deep short-submillimetre observations can probe SMGs too faint 
for selection in longer submillimetre bands (faint SMGs), whose global properties might dif- 
fer from the bri ght SMG populati on (e.g., lower redshift, warmer dust temperatures, lower 
luminosities; see lKhan et al.ll2007l ). Given the paucity of 350 /zm-selected sources, the most 
efficient way to characterise the nature of the population is through deriving source counts 
and analysing the models that best fit the data. This complements the mul tiwavelength anal- 



Khan 


2006; 


Khan et al. 


2007) 



In this paper we present the first constraints on the source counts at 350 fim. We outline 
our observation programme design, data reduction and analysis methodology. We discuss 
the criteria for selecting candidate 350 fim sources and the determination of the number of 
spurious sources in the map. From this we derive the measured source counts from the 
survey. We discuss how the counts reflect the nature of our sources and conclude with the 
implications for future blank surveys in the short-submillimetre wavebands. 



2. OBSERVATION PROGRAMME 



Submilli metre surveys have followed t hree approaches: using gravitational lensing around 



clusters (e.g. ISmail. Ivison fc Blain 



sourc es (e.g., high redsh ift quasars; 



1997), selecting fi elds surrounding known high redshift 



Ivison et al.l 120001 ). and targeting a region of blank sky 



!e.g.. iHughes et al.lll998l ). For a given integration time, the number of detected sources will 
be higher in a lensing cluster survey as compared to a blank survey due to the brightness 
magnification. However, this approach is highly dependent on the cluster mass distribution, 
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which can produce significant systematic uncertainties on the luminosity function of the 
detected population and its evolution. Even in the best possible case (a smooth cluster), 
imprecision in the cluster model could still dominate over the behaviour of the source counts. 
Sub millimetre surveys centred on known high redshift sources run a risk of being redshift 
biased, since these are typically found at other wavelengths, and many are lensed. Addition- 
ally, correlation analyses show a high er probability of fin ding; enhanced source counts over 



typical survey sizes in such areas (see lLagache et al.ll2005l and references therein). 



To avoid the uncertainties associated with the biases listed above, we have chosen to 
pursue a blank survey, which can be implemented in ways ranging from deep, small area to 
shallow, large area surveys. In order to maximise the number of detections in the survey, we 
could discriminate between the two approaches using the following argument: the differential 
number versus flux relationship can be approximated locally as 

N(S) » k d (S/S )-^ [sources sr" 1 mJy" 1 ] (1) 

where N(S) describes the overall surface density of galaxies as a function of flux density S. 
For a given limiting flux, S m i n , the number of sources is 

f°° dN 

N{> S min ) = / — dS N{> S mm ) ex S l ~Z- (2) 

During a single pointed observation, the noise is expected to integrate down as 1/ y/t. Hence 
the number of detected sources, N, is related to the integration time t via 

iVdeep OC t^ 2 . (3) 

If the integration time was instead sub-divided into an equal number of shallower observa- 
tions, this would yield 

N wide oc t. (4) 

Therefore a deep pointing yields more detections per exposure compared to a wider, 
shallow survey, as long as the flux density sensitivity remains at a level where 7 >3. 
For a non-evolving Euclidean universe 7=2.5, but current constraints on the submillime- 
tre g alaxy population sh ow evolution (7 >2.5) for a broad range of brighter flux densities 



[e.g.. iCoppin et al.ll2006l ). Constraining 7 through direct observation requires the detection 
of tens of sources at 350 jttm - a huge demand on telescope time. 

Rather than parameterise the source counts from very small data sets, a more practical 
approach is to discriminate between existing source counts models, in particular those that 
successfully reproduce the IR-submillimetre counts. Using the models in the literature at 
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the time of the sur vey (jFranceschini et al.l 1 1994 ; iGuiderdoni et al.l Il998bl ; iPearsonl 12001 



Takeuchi et al.ll200ll ) the target la sensitivity was based on where the models begin to show 



significant deviations in their source count predictions, with the majority of models having 
7 > 3. This threshold was la=5mJy. 



2.1. Observations 

SHARC II is a 350 /im-optimized camera built around a 12 x 32 element close-packed 
bolometer array. It achieves a point-source sensitivity of ~ 1 Jy s 1//2 in good weather. The 
384 pixels of the SHARC II array image a region of around lfO x 2'.Q on the sky. Its filled 
absorber array provides instantaneous imaging of the entire field of view, sampled at roughly 
2.5 pixels per nominal beam area. 

The 350 /im window is a difficult one for observers: the in-band atmospheric opacity r 
is rarely < 0.8, with signal-to-noise S/N oc e~ T /v / l — e~ T , making efficient observations ex- 
tremely weather dependent. For ground-based far-IR/submillimetre observations, the vari- 
ation in atmospheric emission is the dominant noise source over all temporal frequencies. 
Although rapid image differencing, commonly called chopping, is used to remove the atmo- 
sphere signal, this technique can give rise to a V2 increase in noise and a loss of observing 
time from a chopping duty cycle of < 1. Furthermore, chopping does not adequately remove 
portions of the atmospheric signal that vary fas ter than the chop frequency, something that 



our data reduction analysis has shown to exist (jKhanll2006l ). 



The design of SHARC II eliminates the need to chop. Atmospheric noise is spatially 
correlated, implying that the spatial variation in the atmosphere occurs in the line of sight 
of several pixels. By scanning the detector array over the target region, the celestial signal 

- spatially fixed and constant in time - will be mapped by several detector pixels. This 
scanning technique allows the determination of the individual pixel gains and offsets, and the 
removal of the atmospheric signal on all timescales; least squares fitting can also model other 
instrumental contributions, alongside the simultaneous derivation of the celestial sky map 
and associated uncertainty. Although this modelling will induce some covariance between 
adjacent map pixels, this is small compared to the dominant contribution from photon noise. 
As part of the commissioning phase of SHARC II, we tested a number of Lissajous scan 
patterns, typically using smaller amplitude sweeps of about 15 arcseconds in the x-direction 

- perpendicular to the 32 rows - and 10-20 arcseconds in the y-axia^|. This ensured that the 
entire area was well-covered, with substantial redundancy between detector pixels and map 



2 The amplitude-period ratio should not be much larger than 1.4 arcseconds per second. 
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pixels. 



The survey was awarded 12 half-nights of observing time, commencing in January 2003. 
From that, just under seven hours of good quality data (from observations done in reasonable 
weather in January and March 2004) were obtained, centred on the Bootes Deep Field 
Jde Vries et all booj l at position RA= 14 h 32 m 5!75, Dec= 34°16'47"5 (J2000), during the 
runs in January and March 2004. For these data the in-band zenith atmospheric opacity 
( r 350/im) ranged from 1.0 to 1.4, corresponding to a zenith transmission of around 30 per 
cent. The beam profile was measured on known compact sources, and was verified to be 
within 3% of the diffraction-limited bea mwidth of 8 . 5". A ll observations were taken using 
the Dish Surface Optimisation System (ILeong et al.ll2006l ). which corrects for the primary 
mirror deformation as a function of zenith angle, to improve the telescope efficiency and the 
pointing. 



2.2. Data Reduction and Source Extraction 



The data were reduced using the standard CSO reduction software, CRUSH flKovacs 



20061 ) version 1.40a8, using the advised reduction parameters for deep observations. This 
software implements a self-consistent least-squares algorithm to solve for the celestial emis- 
sion, taking into account instrumental and atmospheric contributions to the signal. Forty 
individual scans, each representing approximately ten minutes of integration time and all 
centred on the Bootes Deep Field position, were reduced simultaneously through CRUSH. 
The output, the CRUSH- reduced skymap was calibrated with the flux density and point 
spread function based on observations of Callisto taken throughout the observing period at 
similar elevations (usually every hour). The flux density of Callisto was derived from the 
CSO SHARC II calibrator catalogue. A thorough treatment of the redu ction metho dology, 



with detailed explanations of the reduction parameters, can be found in (jKhanll2006l ) 



For each pixel in the CRUSH-reduced skymap, a least squares fit for a point source was 
determined. From the CRUSH celestial map, for each skymap pixel j, a sub-map comprising 
all pixels within 16.2 arcseconds (or 10 CRUSH skymap pixels), was extracted. The size of the 
sub-map was chosen to provide a good determination of the source and background, but not 
so large as to require a more complicated background model, whereby four parameters were 
fit simultaneously: source intensity, mean background, and both a horizontal and vertical 
linear gradient. The Callisto point spread function (PSF) was then applied to this model in 
a weighted least squares fit — this is roughly equivalent to smoothing the celestial map with 
the PSF (see Figured]). For each pixel, this fit produces an intensity Sj and an associated 
statistical uncertainty, crj, in units of flux density per beam. These values allowed an estimate 
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of the approximate signal-to- noise (S/N), which we refer to as the "significance" using 

6 = § (5) 

This fitting reproduces the known 350 jttm flux densities of standard calibration sources 
to within the calibration uncertainties, but for faint sources, the map noise is the dominant 
uncertainty. 



2.3. Re weighting the map 

In a map with few detections, the expected distribution of £ will be Gaussian, with a 
variance of one, centered on zero. We define 



(6) 

as the rms variation in £. 

For the Bootes data, H = 1.51 - implying further noise terms not accounted by in the 
CRUSH analysis. While it is possible E > 1 could be due to real structure in the maps 
(such as confusion noise - the statistical variation from unresolved sources), this is unlikely 
given the expected number of detections based on the survey sensitivity (using the models 
in Section [2]). Other models to derive an appropriate scaling factor were considered, from 
a simple constant offset to treating the excess noise as additional variance that is added in 
quadrature to the statisti cal uncertainty from the detector noise using maximum likelihood 
statistics (see iKhanl 120061 ) . but an adequate solution was to simply scale the map by S: 

a] = Eaj (7) 

giving a corrected significance of = The magnitude of H appears stable with 

the integration time: real structure in the sky should be \fi more significant for longer 
integrations. For source counts, the systematics associated with this excess noise are small 
compared to Poisson statistics. From this point, £ and a refer to the adjusted values, £' and 
a'. The adjusted significance distribution in the map is shown in Figure [2, alongside the 
corresponding the survey coverage for the adjusted noise (Figure [3]). It is this adjusted noise 
that is used for source extraction^. 



3 Our analysis shows that this reweighting is still required in maps reduced with newer versions of CRUSH. 
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2.4. Extracted Source Properties 



The corrected significance was used to select candidate detections, where |£| > 3. There 
were th ree positive sources that met t he detection criteria, including the previously reported 
SSG 1 (IKhan et al.ll2005l ; lKhanll2006l ). summarized in Tabled] (note: a is scaled by H), and 
two negative. 

The variation of x 2 with source position gives the position confidence contour, as given in 
Table [IJ quoting 3a positional uncertainties (the best-fitting x 2 position will not necessarily 
match the peak S/N position, as illustrated by SSG 3 in Table [p. 



3. Constraints on the 350 ^m-selected population 

The relation between the measured density of sources and the corresponding flux den- 
sities (the source counts) constrains theoretical models of the source luminosity function 
and its evolution. A thorough treatment of the measured counts would include a variety 
of statistical processes (e.g., confusion noise, errors in the map). But the small number of 
detections in this survey means Poisson noise is dominant. 

Even in the absence of real sources (fis), there will be still a statistical chance of detecting 
a source above the £ > 3a threshold. The mean number of these detections in the entire 
survey is called the accidental rate, (also referred to as the spurious source detection 
rate). If the expected number of 350 [xm sources, both real and spurious, is small, then the 
two types of detections can be considered as independent detection processes, giving the 
total number of detections as /i = \i& + Us- 



3.1. Empirical estimate of the accidental rate 

A standard approach to determining the accidental rate is through using the pixel-pixel 
covariance to produce a model for the expected number of connected regions that lie above 
the detection threshold (3a), assuming these covariances are well-characterised. If the map 
noise obeyed Gaussian statistics, the probability of a pixel having S/N > 3a per beam would 
be 0.00135. The approximate number of map beams is 310 (using the Callisto PSF). Thus 
the expected number of accidental sources would be ha ~ 0.4. In the real CRUSH-reduced 
map, however, the difficulty in characterising the noise (Section 12.31) shows that it is not 
Gaussian, which forces use of an alternat ive method for deter mining; ha- One way is an 
empirical approach, similar to that used in ISerjeant et al. J2003h . based on the fact that sky 
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noise is not correlated with celestial position (a, 5) but real astronomical sources are. 

For each raw data scan, a random rotation^ angle was assigned, and the entire dataset 
with rotation angles was passed to CRUSH for reduction. This has the effect of smearing 
the true astronomical sources while keeping the spatially correlated noise intact. The source 
extraction method of Section 12.21 was used to determine the number of candidate sources 
in the rotated maps. In total, 634 rotated maps were generated this way. Although the 
corrupted-astrometry maps have slightly different area-sensitivity coverage than the original 
map, the uncorrupted map is a random sample from this wider ensemble. The original map 
coverage is typical of the corrupted sample. The excess noise E of the original map is also 
within the range found for the corrupted maps (1.23-1.59). 

The corrupted-astrometry maps produce the greatest density of spurious sources in the 
low-coverage, high-noise regions. However, all the candidate 3a sources in Table [T] are in 
the central region, where a < lOmJy). In this region, the spurious source detection rate is 
Poisson distributed with an expectation of 1.09 ± O.OcI 

With three candidate point sources and an accidental rate of 1.09, the true detection rate 
is poorly determined. However, observations at other wavelengths can assist in determining 
which sources are real. Although this introduces a selection bias, it will be small compared 
to Poisson statistics. Two of the candidate sources in Table [1] are 5a detections at 24 /im 
and 20 cm. The probability of accidental detections at 24 /j,m is 0.3 and 3% for SSG 1 and 
2 respectively. At 20 cm the accidental detection probability is 1% for both sources. Given 
these high-likelihood identifications it is unlikely that either of these two are spurious. 

SSG 3 is more proble matic: the sensitivit y of the 24 jum data suffices to detect 850 //m- 



selected galaxies (see, e.g., lEgami et al.l 12004 ). The non-detection of this sou rce at 24 (im 



and 20 cm suggests it is an atypical SMG, possibly at high redshift (see, e.g., Ilvison et al. 



20021 ; Ilvison et al.l 120071 . although without the radio/24 /im identification no photometric 
redshift estimate can be obtained). But with the expectation of 1.09 spurious sources and 
the multiwavelength identifications of SSG 1 and SSG 2, we assume that SSG 3 is least likely 
to be genuine and so exclude it from further analysis. 



4 Rotation angle is a parameter intended to represent the position angle of the SHARC II array on the 
sky. For present purposes, introducing a random value is nothing more than a simple method of offsetting 
the array astrometry from its true value. 

5 The uncertainty is in the measurement of the accidental rate, not the range on the number of accidental 
sources. 
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3.2. Survey completeness 

To determine the survey completeness the two real sources, SSG 1 and SSG 2 (or SMM 
J143206.65+341613.4 and SMM J143206. 11+341648.4), were removed from the CRUSH- 
reduced skymap and a source of random intensity was inserted into the no-source skymap, 
randomly placed over the entire skymap area, A. The simulated-source map was then fit 
as in Section I2.2[ and the fraction of simulated sources recovered at > 3a was determined 
through a Monte Carlo simulation (with the noise scaled by the same S as the original map). 
The completeness against simulated source flux density is shown in Figure [3], for the deepest 
part of the map (er < lOmJy). 



3.3. Source Counts 



The number of sources detected by a survey in area A to depth S > SVmn will be 

roc 

N det = A N(S) x C(S) dS (8) 

where C(S) is the completeness within the survey area. Typical source count models (e.g., 
those given in Section |2J) are well represented by power laws in flux density, as given by 
Equation 1. Setting Ndet = 2, substituting Equation 1 for N(S), and normalising the differ- 
ential counts at So = 20mJy gives kd ~ 0.035 sources arcmin -2 mJy -1 . The normalisation 
at 20mJy gives the least dependence of on 7 for the present survey, for less than 10% 
variation for 2.5 < 7 < 4.0. 

The uncertainties on kd are set by Poisson statistics. For an ob served count o f two 
objects, the true counts are between 0.53 and 5.32 with 90% confidence ( ]Gehrelslll986l ). The 
uncertainty on kd scales directly with these values. Equation 2 allows direct comparison with 
integral count models. We choose 5 m j n = 13 mJy, again minimizing the dependence on 7 for 
the actual survey, and find 0.84l '£[ sources arcmin -2 with S > 13mJy (as shown in Figure 
HJ), quoting the 90% confidence uncertainty. The variation is <5% for 2.5 < 7 < 4.0. 

In a map with few 3a det ections, a careful consideration of the Eddington bias must be 



applied (e.g.. lEddingtonlll913l ). Because there are usually more sources immediately below 



the flux limit than immediately above it, more sources are scattered above this limit, by 
positive noise fluctuations, than are scattered downwards to below it. Therefore, sources 
close to but above the flux limit have measured flux densities biased high, on average. But if 
we assume a form for the source counts, the effect of Edding ton bias is implicit ly corrected. 
However the deboosted individual flux densities are given in iKhan et al.l (120071 ). 
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DISCUSSION 



The derived integral counts are presented alongside a variety of source count models from 
the literature in Figur e HI The models represent two approaches to source count modelling 
- backward evolution ( 



2001 



Pearson! 120071 . IVaccari et al.l 120071. lLagache et al.ll2005l. iPearson 



Rowan- Robinson! |200ll ) and semi-analytic (IGuiderdoni et al.l Il998bl ) (see lHauser fc Dwek 



2001 



for explanation and detailed descriptions of these methodologies). The 350 /im pop- 
ulation, like other submillimetre-selected populations, is evolving with redshift, with num- 
bers more than an order of magnitude higher than no-evolution predictions. At 90% confi- 
dence we are able to reje ct the No Evolution model, as well as the no-ULIRG model from 
Guiderdoni et al.l (Il998bl ). But due to the small sample size the bulk of the 350 /im models 



cannot yet be discriminated or rejected. 

The small area of this survey means the source counts will inevitably be affected by 
cosmic variance. But the number of 5a 2 4 mm d etections within the SHARC-Bootes area 
compared to the counts of IPapovich et al.l (120041 ) suggest an underdensity in this field (see 
Khan et al.l 120071). Also , the photometric redshifts of the two detected sources (z ~ 1 and 
z ~ 2; iKhan et al.ll2007l ) make it unlikely these objects are related to each other. 



For comparison we plot the 450 /im counts from lSmail et al.l (120021 ) in Figure H] assuming 
an Arp220 SED template to transform the 450 /im counts to 350 /im (the 450 /im flux density 
of lOmJy being roughly equivalent to a 350 /im flux density of ~16mJy). Although this is 
a crude shift it appears consistent with the 350 /mi counts. These counts are also consistent 
with t he 350 /tm limits (at ~25 mJy) on 850 /mi-selected sources presented in ICoppin et al. 
J2007h . 



Using the relation of Fixsen et al.l (Il998l) . the 350 /im contribution to the cosmic infrared 
background (see, e.g., lLagache et al.l 120051 ) is 0.65 MJy sr _1 . From the source counts we 
estimate resolving around 30% of the 350 /im background at 13 mJy (with the entire 350 /im 
background being resolved at a flux den s ity of ~0.5 /iJy). Although this is rou ghly the double 



the n umber resolved by the lSmail et al.l (120021 ) survey at S 350 =1Q mJy (see also lLagache et al. 



20051 ). the counts are extremely steep in this flux density domain and thus small increases 
in sensitivity result in large resolved fractions. 



Khan et al.l (120071 ) discuss the spectral energy distributions of the two sources detected 
and show that the luminosities are ~ 1O 12 L and dust temperatures are in the range 30-40 K, 
placing them in the region of luminosity-dust temperature space between local I R-luminous 



galaxi es and the colder, more luminous, an d much mor e mass ive SCUBA sources iBlain et al. 
( 120041 ) . This supports the argument of IKhan et al.l (120051 ) that the short-s ubmillimetre 



might sample a warmer SMG population. Indeed the upper limits at 1200 /im (IKhan et al. 
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20071 ) imply that the SHARC II sources may lie below the detection limit of the SCUBA 
instrument at 850 /im. Given this survey resolves a lar ger fraction o f the short-submillimetre 
background compared to the 850 ^urn-bright sample of ISmail et al.l (120021 ). it is possible faint 
SMGs outnumber SCUBA-bright sources (defining a faint SMG as 5*850 ^$5mJy). 

In order to better understand the nature of the short-submillimetre population it will be 
necessary to increase the number of sources, sampling a larger dynamic range in flux density. 
This can be achie ved through follow-up imaging of SM Gs selected at long-submillimetre 
wavelengths (e.g., iKovacs et al.ll2006l ; ICoppin et al.l 120071 ) . or through deep surveys similar 
to this one. But a far more efficient way will be through space-based and balloon-borne 
surveys. 



ESA's H erschel Space Observatory (Herschel, due for launch in ~2008; iPilbrattl 12002 
Harwitll2004l ) will carry out both medium and deep surveys i n the short-submilli metre wave- 
lengths (250, 350 and 500 /zm) with the SPIRE instrument (IGriffin et al.l Il999r) . Similarly, 
the Balloon-borne Large Area Submillimeter Telescope (BLAST; Devlin et al.ll2004l ) will con- 
duct deep, large area surveys at submillimetre wavelengths including 350 /iin. These surveys 
will select large numbers of sources, making it possible to assess the relative contribution of 
bright and faint SMGs to the short-submillimetre background, and determine, through mul- 
tiwavelength analysis, whether the global properties of the short-submillimetre population 
are different from the SCUBA-bright SMG population. 

But the turnover in the 350 /i m differential counts is predicted to occur in the flu x 
density range 5< S350 <20mJy (e.g.. lLagache et al.ll2004 |Pearsonll2007l ; IVaccari et al.ll2007l ). 
which is below the 20 beams per source confusion limit (~21mJy) for the Herschel SPIRE 
wavebands (BLAST will also be confusion-limited at flux densities <25mJy; iPearsonl 120071 ). 
This is a powerful diagnostic to discriminate both evolutionary models and the sources 
dominating the 350 /im background — sub-L* galaxies that will dominate the CIB and the 
volume-averaged star formation rate — hence ultra-deep ground-based 350 /im surveys could 
be the only plausible opportunity to detect this break for the foreseeable future, with the 
same argu ment applying to su rveys in other short-submillimetre bands, e.g., 450 /im with 
SCUBA 2 ^Holland et allbood ). 



5. CONCLUSION 

The SHARC-Bootes survey is a ~6 arcminute 2 blank field survey that achieves a peak 
la 350 /im sensitivity of ~5 mJy Having accounted for artificial sky structure in the map, 
we detect three candidate sources with S/N> 3a. From our three detections, we use a Monte 
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Carlo simulation to deduce a spurious source detection rate, which is Poisson distributed with 
an expectation of 1.09 within the central region of the map. Deep 24 /zm and 20 cm imaging 
is used to confirm the detections and exclude spurious sources. From this identification in 
other bands, and with a likelihood of one source being accidental, we believe there are two 
real 350 /xm-selected sources in our survey. 

Our source count indicates that the IR-luminous population at 350 fim is evolving with 
redshift, with the no-evolution scenario rejected at 90% confidence. 350 fim surveys with 
BLAST, and after that, Herschel, may be unable to probe sources below our current survey 
detection threshold (due to the constraints of source confusion) where the differential counts 
are expected to turn over, therefore future ground-based observations should be designed to 
constrain this break through ultra-deep surveys. 
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Candidate Source 


Flux density [mJy beam 1 ] 


Peak significance [a] 


Position [J2000] 


SSG 1 


23.2±6.5 


3.6 


14:32:06.65±0.26 +34:16:13.4±3.4 


SSG 2 


17.1±5.4 


3.2 


14:32:06.11±0.28 +34:16:48.4±3.2 


SSG 3 


19.9±7.1 


3.0 


14:32:07.46±0.39 +34:17:19.3±8.1 a 



Table 1: Flux densities at the best-fitting position for sources in the SHARC-Bootes Field 
with significance >3.0a (quoting 3<r positional uncertainties). 



"This confidence region is affected by the close proximity of two spurious sources with -3.6 and -2.2a. 




Fig. 1. — Point source intensity top left and noise maps top right from the least square fit to 
the CRUSH reduced map [Jy per beam]. The bottom panels show the significance map left, 
and the map pixels with £ > 2.8 right. 
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-4-3-2-101234 

Significance {a} 

Fig. 2. — Histogram of the uncorrected significance values, £ red line (in units of a), and 
corrected values, £' 6/acA; line, in the SHARC-Bootes map. Scaling by H reduces the rms of 
the unsealed distribution to be ~1. 
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S {mJy} 

Fig. 3. — Dashed-line: adjusted coverage map [arcmin vs la noise in Jy per beam] for 
the SHARC Bootes survey, after a is scaled by S. Thick-line: the completeness for the 
SHARC-Bootes survey, showing the fraction of simulated sources (with truth flux density, 
Stme [ m Jy]) recovered at > 3cr. We only consider the sources with la noise <10 mJy (A deep = 
3.51 arcmin 2 ) thereby excluding the map edges. 
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Fig. 4. — The 350 /im source count fro m this survey with corresponding 90% confidence lim- 
its, and the shifted 450 /im counts from lSmail et alJ ( 2002 ). The horizontal error bars on the 
450 fim counts give the range on the Arp220 SED at 1 < z < 2, with the actual flux density 
being the mean. Thse counts are plotted with a selection of popular source count models 



from the literature , representing two methodologies of count modelling: 



Dackwards evolu- 



tion ((Pearsonl 120071, IVaccari et al.ll2007l. lLagache et al.ll2004 



Pearson 2001 



Rowan- Robinson 



200ll ) and semi-analytical (jGuiderdoni et al.lll998bl ). The lGuiderdoni et al.l (11998bl ) A model 
does not include a ULIRG population, whereas model E includes a strongly evolving ULIRG 
component w ithin the extragalactic population. Also shown is the No Evolution model of 
Pearson 2007. 



